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Abstract. This paper is devoted to presenting in detail the non-perturbative 
renormalization group (NPRG) formalism to investigate out-of-equilibrium 
systems and critical dynamics in statistical physics. The general NPRG framework 
for studying non-equilibrium steady states in stochastic models is expounded and 
fundamental technicalities are stressed, mainly regarding the role of causality 
and of Ito's discretization. We analyze the consequences of Ito's prescription 
in the NPRG framework and eventually provide an adequate regularization to 
encode them automatically. Besides, we show how to build a supersymmetric 
NPRG formalism with emphasis on time-reversal symmetric problems, whose 
supersymmetric structure allows for a particularly simple implementation of 
NPRG in which causality issues are transparent. We illustrate the two approaches 
on the example of Model A within the derivative expansion approximation at 
order two, and cheek that they yield identical results. We stress, though, that the 
framework presented here also applies to genuinely out-of-equilibrium problems. 
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1. Introduction 

The theoretical understanding of out-of-equihbrium systems, and in particular 
non-equilibrium steady states, is nowadays one of the major quests of statistical 
physics. These systems exhibit a great variety of behavior, including genuinely 
non-equilibrium phase transitions and critical phenomena, ubiquitous strong-coupling 
regimes, recurrent absence of lower critical dimension, etc. Numerous partial results 
have been obtained at various levels, but all-purpose, analytical, and systematic 
methods to treat non-equilibrium problems are scarce. 

The nonperturbative renormalization group (NPRG) is one such method, and this 
paper is devoted to detail its implementation for out-of-equilibrium systems. Indeed, 
many remarkable results have already been obtained in statistical physics using the 
NPRG approach not only at equilibrium in systems such as frustrated magnets [1], 
the random field or random bond Ising model [2] , membranes [■]] , bosonic systems [4] , 
but also out-of-equilibrium, where one can mention important advances in reaction- 
diffusion systems [5], and for the Kardar-Parisi-Zhang universality class [()]. 

The starting point for applying NPRG methods is a field theory. Such a field 
theory is usually obtained by following one of two possible routes. Either one starts 
from an effective mesoscopic description in terms of a Langevin equation and one then 
resorts to the Janssen-de Dominicis formalism [7] to construct a functional integral 
representation of this equation, upon introducing an additional Martin-Siggia-Rose 
response field [8]. Or, and this is the other route, one starts from a microscopic 
master equation that can be mapped onto a functional integral following the Doi- 
Peliti procedure [9] . When both descriptions coexist the resulting field theories arc of 
course equivalent (see however [10] for important subtleties in the derivation of the 
equivalence of the Janssen-de Dominicis and Doi-Pcliti field theories) . 

Both field theoretical formulations are derived in discrete time, and a crucial 
ambiguity arises when performing the continuous-time limit [11, 12]. This problem 
implies the choice of a prescription in the continuum version (usually related to the 
Ito's or Stratonovitch's discretization scheme), the consequences of which still deserve 
discussions {e.g. [13]). Roughly speaking, the Janssen-de Dominicis transformation 
for instance involves a determinant which depends on the original discretization 
of time. This determinant can be conveniently set to unity upon adopting Ito's 
discretization scheme and imposing the related prescription in subsequent calculations. 
This prescription is easily implemented in perturbation theory: it amounts to cancel 
out tadpole diagrams [12]. 

The implications of this choice of time-discretization within the NPRG formalism 
have yet to be clarified and is one of the main subjects of this paper. We shall first 
identify the consequences of Ito's discretization in NPRG flow equations and then 
device a simple procedure to systematically encode them. This will be the opportunity 
to review some other important technical aspects arising when implementing, out-of- 
equilibrium, the most popular approximation scheme of the exact NPRG equation, 
the derivative expansion [14, 15]. 

To guide us in the task of understanding the consequences of Ito's discretization 
in the NPRG framework, we first review, in section 2, the discrete-time version of 
non-equilibrium field theories ensuing from the two standard Janssen-de Dominicis 
and Doi-Peliti procedures, analyze causality issues, and recall the usual prescription to 
encode Ito's discretization in perturbation theory. Section 3 is dedicated to expounding 
the specificities of the NPRG framework when applied to non-equilibrium systems. 
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Causality issues arc postponed to Section 4, which focusses on identifying the relevant 
prescription associated with Ito's time-discretization in the NPRG context, and where 
we propose a simple "regularization" procedure to encode it. Even though it can be 
(and has been) applied to genuinely out-of-equilibrium problems, we illustrate it, in 
Section 5, to a relaxation-towards-equilibrium problem, the critical dynamics of Model 
A. In Section 6, we show that for such models satisfying time-reversal symmetry, one 
can build a supersymmctric formulation of NPRG. Not only this formulation is elegant 
and convenient for actual calculations but, more importantly here, it is free from any 
ambiguity regarding the continuous-time limit [12]. Once the supersymmetric NPRG 
formalism is set up, it is applied to Model A, which allows us to check that the two 
approaches yield identical results. 

2. Field theories for nonequilibrium systems 

As it constitutes the basis of the following discussions, we start by recalling the 
derivation of the field theory associated with an out-of-cquilibrium stochastic model. 
We review the two standard procedures, starting from a Langevin equation or starting 
from a master equation. In both cases, the choice made for the continuous-time limit 
plays an important technical role both perturbatively and nonperturbatively. 

2.1. Doi-Peliti functional 

Many models (such as reaction-diffusion systems) arc defined at a microscopic scale 
by a set of dynamical rules. The time evolution of these systems is governed by a 
master equation, which, in an operator representation, can be written as 

dt\P{t))^W\P{t)) (1) 

where \P{t)) ~ J2n^i^^^)\^) associated with the probability law of the 

system, \n) are the Fock states (number states) and VV = yV{a^,a) is the operator, 
expressed in terms of creation and annihilation operators, representing the kernel of 
the master equation. Eq. (1) can be cast into a field theory following the Doi-Peliti 
procedure ['■)]. We do not review it in detail here (see e.g. [9, 16]), but merely stress 
a few points regarding the ordering of discrete times. 

The mean value of an observable O at time t can be expressed in the operator 
representation as 

{0){t) = {V\de^'\P{Q)) (2) 

where (P| is the projection state satisfying {V\P{t)) = 1 and (T'lW = 0. The 
procedure consists in splitting the evolution operator exp(yVt) in + 1 time slices of 
width T and in inserting in each a closure relation written in terms of the coherent 
states 

{0){t) =j(jl ^) {V\d\M{Me'^l<I^N-i) ■ ■ ■ (0i|e^n0o)(0o|P(O)) (3) 
where (P4'a = dTm{4>a)dTZe{(j)o.). At first order in r and after re-exponentiation, it 
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writes 

By computing the usual nonvanishing overlap {(j}k\4'k-i) of coherent states, one 
deduces (at first order in r) the generating functional 

•' Q=0 V k=l fc=l / 

(5) 

where the initial conditions and final time contributions have been omitted since our 
focus is on stationary states and a loss of memory of initial conditions is assumed (see 
e.g. for details). Note that the two fields (pk and (pl are complex conjugate and 
that we have chosen real sources jk^jk (see Appendix 8.1). 

The continuous-time limit consists in taking the limit r — > 0. In this limit, the 
discrete sums become integrals over time and 4'ki't'k ~ 4>k-i)/T tends to (jfOtcj). As 
we shall see in the following, this continuous version of the model becomes ambiguous 
when the inverse of the operator dt5{t — t'), that is 0{t — t'), appears at coinciding 
times. To remove these ambiguities, we need to refer to the original time-discretized 
version of the model, Eq. (5). Let us already mention that the way the coherent states 
were introduced in the derivation above leads to the time ordering of Eq. (5) where 
the 0* field always appears at a time equal or larger than the time of the ((> field. We 
show in the following that this also occurs for field theories derived from Langevin 
equations when Ito's time-discretization is chosen. 

An example of a field theory derived from a master equation such as Eq. (1) is 
that of branching and annihilating random walks in the universality class of directed 
percolation [ Hi]. A particular model belonging to this class consists in a set of identical 
particles A diffusing on a lattice at rate D and undergoing, upon encounter, the 
reactions: A — > 2A and 2A — > with rates a and A and, possibly, A — > with rate 
jji. These models exhibit a phase transition when the system falls into the absorbing 
(empty) state which sets in, e.g., when the branching rate cr is decreased. Obviously, 
this dynamics is irreversible since detailed balance is violated. The transition is thus 
genuinely out-of-equilibrium. 

The reaction-diffusion system above can be studied using the Doi-Peliti formalism. 
Its action writes (in the continuous-space and time limits and after a shift and a 
rescaling of the fields and (/>*): 

= j {dtc^ - DV''cf>) + [ii - a)<i>* <i> + V2^(j)*<i>{(j) ~ (j)*) + K<t>* (1^?) ■ 

(6) 

All field theories derived from this formalism share some common features: (i) 
their actions are proportional to the field 0* (which is an important property, see 
Section 4.2), (ii) they are very similar to the corresponding field theories derived from 
Langevin equations (see [10] for a discussion about the correspondence between the two 
formalisms) . We now review the derivation of field theories ensueing from Langevin 
equations. 
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2.2. Janssen-de Dominicis functional 

Suppose that a field (j){t, x) evolves according to a stochastic dynamics described at a 
mesoscopic level by the Langevin equation: 

dtcj){t, x) = -niF{dp{t, x)) + N{dp{t, x))Cit, x) (7) 

where F represents the deterministic part of the evolution which can depend on (f> and 
its space derivatives, fli denotes a constant and uniform relaxation rate, and <^(t, x) 
is a Gaussian noise: 

P(C) = ^^e-j/*,.-?'/"^ with / =fd'^xdt (8) 

so that 

(C(t, f )C(i', a?')) = 2^(2 5{t - t')5\x - x') . (9) 

One wants in general to compute averages of functions of the field 0{4){t^ x)) over the 
noise distribution: 

{Om^ j VQPiOOic^^) (10) 

where (/'^(i, x) is the solution of Eq. (7) for the realization C, of the noise. This average 
can be conveniently written as 

^jvc P{0 J vcj, 6{dtc^ + n,F{(t>) - 7V(0)C) o{4>) (n) 
= / PC ^(C) / V(j)V[i^] ^m^+i^iFM-NWQ j^^^ 

where the functional analogue of the usual identity 

5{x-xo)=5{f{x))\f'{xo% (12) 

with f{x) assumed to have a unique zero xq, has been used. In Eq. (11), the Jacobian 
writes: 

y,„.|aet(a + n,^M_«)c)|. (13) 

Notice that to go from the first to the second equality in Eq. (11), the initial conditions, 
which arc in principle necessary to solve the Langevin equation, have been omitted. 
As in the Doi-Pcliti formalism, it is implicitly assumed that the dynamics under study 
leads at large time to a steady states with a loss of memory of the initial conditions. 
This derivation is therefore only valid in (or sufficiently close to) the stationary state. 

At this stage, two routes can be followed. The first route consists in introducing 
a set of (conjugate) Grassmann fields rj and f] to exponentiate the determinant as 

J= /"pr^2?r^e-'''.-"(^'+^^^^-^^)'' (14) 



(assuming that the absolute value plays no role which is the case if there is only one 
solution to Eq. (7)). This route is a priori complicated since it doubles the number 
of fields. It is in fact convenient only when the system relaxes to equilibrium. In this 
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case, the model satisfies a supersymmetry that allows one to reformulate the whole 
theory in terms of a unique superfield which makes the corresponding formalism rather 
simple. We deal with this case in Section 6. 

The second route exploits the identity det = expTrln, so that the determinant 
comes as an additional term in the exponential in (11). This term appears to be 
proportional to the inverse of the operator {dt)5{t — i'), which is 9{t — t'), evaluated 
at t = t' [12, 13]. This constant 9{0) = e is ill-defined because it depends on the 
precise (discrete) ordering of times [11]. The simplest choice e = leads to = 1 
and corresponds in discrete time to Ito's discretization choice. Once the Jacobian is 
set to unity, one can finally integrate in (11) over the Gaussian noise distribution and 
deduce the generating functional for correlation and response functions 

Z[j,j] = J V(j,V[i4>] e-'S[^.^l+/*,^J>+J^ (15) 

with 

5[0,0]= [ 4>{dt(t> + niF{4>))~^2N^{4>)4'^. (16) 

Jt,x 

As in the case of reaction-difl^usion systems, and for the same reasons, the 
ambiguity on the value of 0(0) has to be removed both perturbatively and 
nonperturbatively to make consistent calculations from Eqs. (15,16). Since this 
ambiguity is related to the continuous-time limit we first review the discrete time 
version of the action (16). 

2.3. Field theory from discrete-time Langevin equations 

The Langevin equation (7) in Ito's discretization reads (assuming for simplicity 

^1=^2 = 1) 

= rF(0„_i) + 7V(</)„_i)C„-i (17) 

where n, n — 1 are time indices and where r denotes the time step. Any reference 
to the space dependence has been omitted here since it plays no role in the following 
discussion. With this stochastic process is associated a transition probability 

T^((/)„|(/)„_i) = (0„-i)t) exp I ^2IP^ — i> I ^ ' 

that can be rewritten as the following Gaussian integral 



JV((^„|0„-i) = / -^--^T{<i>n-,<\)n\<\>n-\) (19) 

with 

^Mn, = exp (-^„(0n - K-\ - Ti^('/>n-l)) + riV^ (0„- 1 )0,^j) (20) 

where </>„ is the conjugate, or response, variable associated with the transition from 
</>„_! to (/)„. For a Markov chain of N transitions between times and tj<! = tN , the 
total transition probability is the product 

JV 

Pn{^\M\ • ■ • ; 0JV,ijv|'/'O, 0) = T.r(0ri|0n-l)- (21) 
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Hence, one obtains the generating functional of correlation and response functions 
[11]: 



N 

Zm J ■'-J- 



1^ 



n=l 

N 



with S[cl), 0] = T 0„ ((0„ - 0„_i)/t - F(0„_i)) ~ iV^(0„_i)0^. 

n=l 

(22) 

We emphasize that in this action all the fields appear at a time larger or equal 
than the times of all the (j> fields as in the Doi-Peliti procedure for reaction-diffusion 
systems, Eq. (5). The ambiguity in the value of 9{0) is therefore the same in both 
cases. Before embarking into the NPRG formalism let us first recall (and reformulate 
in a way useful for NPRG) how Ito's prescription is dealt with in perturbation theory. 

2.4- Dealing with Ito's discretization in perturbation theory 

To identify the relevant prescriptions, we use the time-discretized version of the field 
theory (22), which is the only form free from ambiguity. As mentioned above, the 
time ordering and the structure of the Janssen-de Dominicis and of the Doi-Peliti 
functionals are the same. Hence the considerations developed below hold for both 
functionals. 

We consider the quadratic part of the action 

N 

So[<j), 0] = r ^ 0„ {4>n - (pn-l) It + W(t>n4'n-1, (23) 
n=l 

with typically w ^ V'^ + m^, and denote Zq the corresponding generating functional. 
Integrating over in Zq produces a product of delta functions which enforces the 
relation 

M}{^o...^N) = 'Cjo.-.jN) (24) 

where *(•..) is the transpose operator, and M is the N x N matrix with components 
Mij = Sij + wSi-ij. This matrix M is precisely the matrix whose determinant 
appears in Eq. (13). It is clear from its definition that dct(M) = 1 and, since 
this result straightforwardly generalizes to any interacting theory, this proves that 
in Ito's discretization the Jacobian is unity. Inverting this matrix: M~j^ = Sij + 

"^^^i {—w)''Si-k,j allows one to perform the remaining integrations over the and 
one finds 

Zo[{jnJn}] = exp I M'l jn j . (25) 

\m,n— 1 / 

One can then compute the free two-point response function (the bare propagator Gq) 

Go,ki = (Ml) = = Mki' = 0{l < k){~w)''-' (26) 

which vanishes if the field (f> appears at an earlier time than the field </> as a signature 
of causality. It is important to notice that {4>k4'k) = 1 whereas {(j)k4>k+i) ~ which is 
the reason why the continuous-time limit t — >■ is delicate. 
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Figure 1. Example of a pcrturbativc tadpole diagram. The time grows from 
left to right and the arrows go from ij> to 0, allowing to keep track of causality. In 
this example, the interaction term in the action is ifi^clP and the time separation 
between the and the </> fields at the vertex is the time step t of Ito's discretization. 
The graph is vanishing because the arrow in the tadpole is going backward in time. 



In perturbation theory, the ambiguity occurs only in tadpoles in which the bare 
propagator Go(t — t') oc 0{t — t') starts and ends at the same vertex. The value of 
a tadpole is therefore proportional to G(0) oc 9{0) which is ill-defined. Bearing in 
mind that in discrete time the interaction vertices come with the fields shifted in time 
according to Eq. (22), that is, are of the form 4>^(j)'^_i, one concludes that the bare 
vertices actually involve two times separated by the time step r. As a consequence, 
the propagator appearing in a tadpole - joining a (f> with a - is thus {<j)k-i<j)k) , 
which is vanishing because of (26), and not {4>k<t>k) (see Fig. 1). This result, grounded 
on the analysis of the discrete time field theories, implies that in the continuum the 
propagator in a tadpole is ((/)(t)0(t+T))o = Go{t—{t+T)) = Go(— r) and that therefore 
9{0) should be understood as: 

9(0) = lim 9{-t) = 0. (27) 

This yields the standard prescription used in perturbation theory: in Ito's 
discretization, the 9(0) appearing in Go(0) must be set to in tadpoles. 

Let us notice that the need for a prescription in perturbation theory is directly 
related to the fact that in continuous time the bare vertices arc local in time and thus 
are singular functions of their time arguments (Dirac functions such as S{t — t')). The 
interesting point for what follows is that, since only the product of a bare propagator 
and of a bare vertex matters, one can remove the ambiguity, either by considering 
that in the vertices the (f) fields are shifted in time with respect to the fields, or 
by considering that this shift occurs in the propagator. In this second solution, Ito's 
prescription Eq. (27) is effectively implemented by shifting the time t' of (j){t') in 
Go{t~ t'), that is, by replacing {(f){t',x')(j){t,x))o by 

(0(t', f )0(t, x))o,e = m' + ^')^it, ^))o (28) 

with e — > 0^ or, equivalently, in Fourier space 

{4>iuj', q'^u:, q))o.e ^ e"""' (^(c.', g')</.(c., q))o. (29) 

Doing so precisely amounts to changing in tadpoles Go(0) by Go(— e), that is, 9{0) by 
9{~-e) as it should according to Eq. (27). 
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3. The non-perturbative renormalization group for out-of-equilibrium 
models 



The general idea underlying the NPRG is the same for equilibrium and out-of- 
equilibrium problems: one builds a sequence of effective models that interpolate 
smoothly between the micro- and the macro-physics and that consist, at an 
intermediate length-scale, in the integration up to this scale over the stochastic 
fluctuations. Two important features of out-of-equilibrium problems must nevertheless 
be specifically considered, the doubling of each field by a response field, and causality. 



3.1. General NPRG formalism for out-of-equilibrium models 

We consider the field theory given by Eq. (22) in Ito's discretization (or Eq. (5) for 
reaction-diffusion systems). As we shall show in the following, this time-discretization 
requires to implement consistently a prescription on the propagator of the theory. We 
ignore this for the moment as we will show in Section 4 that it can be easily taken 
into account a posteriori. 

As in equilibrium, one wants to generate, through progressively averaging 
over fluctuations, a sequence of scale-dependent models, whose effective actions r„ 
interpolate smoothly between the bare action S of the initial model and its effective 
action F, that is, the (non-equilibrium analogue of the) Gibbs free energy. This 
requires that at the microscopic scale A (the inverse of the lattice spacing for instance) 
Fk^a = S and that at scale k — , Fk=o = T. At an intermediate scale A > k > 0, 
the slow modes (with respect to k) are decoupled by construction in order to make F^ 
(almost) equal to the effective action of the rapid modes. 

The decoupling of the slow modes is achieved by adding a quadratic (mass-like) 
momentum-dependent term ASk, to the original action. The required properties of 
this term are (i) to be very large at scale k = A for all momentum modes (so that 
all stochastic fluctuations are frozen and the mean field becomes exact: T\ = S - see 
Appendix 8.1), (ii) to vanish at scale k = (so that the original model is recovered: 
Fk=o = r) and (iii) to give a "square mass" of order to the slow modes when 
A > K > in order to decouple them from the long distance physics (integration 
over the rapid modes only). We therefore build, as in equilibrium, a scale-dependent 
generating functional of the correlation and response functions: 

Z.[jJ] = Jvcl,V4>cxp (^-S -AS^ + J * J(x).$(x)^ (30) 

with X = (i, x) 

*(x, . ( ) ,/(x, = ( ^ ) ,31) 

and 

= *$(x).i?4x-x').<i>(x') (32) 

where Rk is the 2x2 matrix of mass-like cutoff functions that suppress the fluctuations 
of the slow modes (see Appendix 8.2 for the definitions used throughout the text). 

We also define as in equilibrium the generating functional of connected functions 
W^si'^] ~ log^K[J] and its (modified) Legendre transform F„ by 

r«[*]+iy«[J] **(x).i?4x-x').*(x') (33) 
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where 

= (34) 

It is convenient to define two notations for vertex (one-particlc-irredueiblc) correlation 
functions: 

and 

4"'")[{x,}, {x' ■}; *] = (36) 

5V(xi)...<5^(xVO 

where {x^} stands for (xi, . . .x„) and {x'j} for (x'l, . . . ,x'ft) (see Appendix 8.2 for 
the precise correspondence between both notations). 

The exact flow for is given by Wettcrich's equation [17, 14] 

9Kr«[*] = iTr/ a,i?,(x-x').GK[x,x';*] (37) 

^ Jx.x' 



with Gk — 



" ("21 

the full field-dependent propagator and where is the 

2) 



2x2 matrix whose elements are the F,..,,. 

The flow equation of the two-point functions evaluated in a uniform (in time and 
space) field configuration is also particularly important. It is trivially derived from 
Eq. (37) and can be written as, in Fourier space 

5,Fi%(p;vl/„)=.Tr / a,i?(q).G(q). 

X (38) 
ff ) (p, q).G(p + q).f f ) (-P, p + q) - - vf (p, -p, q) .G(q) 



where p = {^tP) represents a couple (frequency, momentum). The dependence is 
implicit in the right hand side of Eq. (38) to alleviate the notations and f P"* (resp. 
f|^'*) is the 2x2 matrix of functional derivatives of f f^-* w.r.t. '^i (resp. w.r.t. 
and \l/j) evaluated in the configuration ^E'm, (sec Appendix 8.2 for more details and 
notations) . 

Of course, as in equilibrium, such a flow equation is not closed, and is only the 
starting point of an inflnite hierarchy. Approximations are needed to obtain any 
useful result. The most widely used so far is the derivative expansion, where the 
vertex functions arc, in Fourier space, expanded as a power series of their momenta 
and frequencies and which we will use in our application of NPRG to Model A in 
Section 5. 



4. Ito's prescription and NPRG 

The flow of a vertex function involves integrals of products of other vertex functions 
Fk"'""* and of full propagators G^ (see e.g. Eq. (38)). As in perturbation theory, 
an ambiguity can only appear when two times coincide. If the (full) vertex functions 
Fk"'"^ are smooth functions of their time arguments, no ambiguity can arise because 
the times t and t' of the two ends of a propagator Gre(t — t') joining two legs of F^"'"'' 
is integrated over so that the value of G„(0) becomes immaterial in these integrals. 
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Once again, this value plays a role only if the propagator Gi^{t — t') is multiplied by 
a singular function such as, for instance, a Dirac function d{t — t'). 

This is precisely what happens within the derivative expansion where the vertex 
functions are, in direct space, expanded as (formal) power series of Dirac functions 
and of their derivatives: (5'™^(t — t'). As in the perturbative case and for the same 
reason, the trick to get rid of the ambiguities is to shift in the full propagator Gk. 
(actually in Gk^'^'' = (0^) which is the only causal two-point function) the time of the 
(j) field by an infinitesimal amount e. Before proceeding, let us prove a theorem. 

Causality of the response functions and choice of cutoff functions 

The causality property of Eq. (26) can be straightforwardly generalized to all response 
functions both for the free field theory and for interacting theories truncated at the 
mean field (tree) level [11]. In discrete time it writes: 

n -^fe. Ilk)^^ if 3 Ij/ I, > h (39) 

for n > and n > 0. This equality means that for a response function to be nonzero 
the largest time must be that of a field. We now prove (in continuous time) that 
this equality holds at any scale k of the NPRG flow. In fact, to preserve this property 
requires to appropriately choose the cutoff functions i?^ and we now discuss this point. 

Most often, it turns out to be sufficient for the decoupling of the slow modes to 
choose an anti-diagonal matrix i?^, that is, a A5k term that couples only (p with (p. 
The exception to this rule arises when the symmetries of the model, that should of 
course be preserved by the cutoff term, enforce the presence of an additional (diagonal) 
term, proportional to (j)^ (as for instance occurs in the field theory associated with the 
Kardar-Parisi-Zhang equation [18, 6]). To the best of our knowledge, a cutoff term 
proportional to (j)^ is never necessary and we do not consider it in the following (it 
would become problematic for the proof of the property (39) and of the fact that 
is proportional to i), see Section 4.2). 

The cutoff term coupling with (j) must of course also preserve causality. 
It is therefore either a function "independent" of time (actually of frequency): 
RK,{t — t',x — x') — )■ 5{t — t')Rii{x — x') or proportional to 0{t — t'). For the sake 
of simplicity we consider in the following a "time-independent" cutoff function but 
the generalization to a general causal cutoff would be straightforward. To further 
simplify, we consider only a AiS^ term that couples 4> with (j) although the following 
arguments can be rather straightforwardly generalized to an additional cutoff term 
proportional to (jP' . We finally choose a cutoff term with a (j) shifted in time: 

A5„= f^R^{x^x')(t){t,x)^{t + e,x'). (40) 

The proof that (39) is preserved at every scale k is made by induction. First, 
the property is obviously satisfied at scale k = A since the mean field is the initial 
condition of the NPRG flow: Fa = 5 (for a proof of this equality, see Appendix 8.1). 
Let us suppose that it holds at a scale k. Then, at scale n—dn the property is preserved 
if the variation of the Green functions coming from the RG flow also satisfies Eq. (39). 
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The starting point of the proof is the NPRG flow equation for that writes: 

d^W^ = [ _ dMS-x') {4>{t,x)4>{t + e,x'))^ (41) 

where the index k in (. . .)„ means that the average is taken in presence of the AS^, 
term. The flows of the connected fmrctions follow from Eq. (41): 

(5"+"w r _ _ 

9k- ^= {(j,i...cbn(l){t,x)^{t + e,x'))^ (42) 

OJl--- OJn Jt,x,x' 

where the indices i stand for {ti,Xi). By hypothesis, the function appearing on the 
right hand side of Eq. (42) is nonvanishing only when its largest time is that oi a. (j) 
field. It cannot be (j){t,x) since 4>{t + e^^?') is posterior and it is therefore one of the 
other (pi fields. This proves that the contribution to W^"'"^ of the momentum shell 
dn is nonvanishing only if its largest time is that of a field. By iteration from the 
inital condition we conclude that the property Eq. (39) holds for any n. 

Note that a similar result can be derived for the one-particle-irreducible vertex 
functions Fk"'"' with the difference that the latest time must be that of a -0 field (this 
subtlety coming from the fact that these functions arc amputated of the propagators 
of their external legs). 

^.2. Ifb's prescription on the full propagator Gl^'^^(t,t') and on 

The previous result implies in particular that the running full (connected) propagator 
G^K'^\t,t') = {4){t' ,x')(t>{t,x)) = Wi^''^\t,t') remains proportional to e{t - t') all 
along the flow since it is nonvanishing only when t is larger than t' . The ambiguity at 
coinciding times remains therefore identical to the one encountered pcrturbativcly. As 
emphasized in Section 2.4. the ambiguity can be equivalently removed in perturbation 
theory by shifting the fields at the vertices or in the bare propagator. This 
second way of shifting time can in fact be formulated at the level of the discrctized 
theory and not only in the perturbativc expansion, therefore it can be exploited 
also nonpcrturbativcly. The way to keep track of a shift of time in the propagator 
Gl^'^^tjt') while Fourier transforming is to modify it as was done in the perturbative 
case on Go{t,t'), Eq. (29). The nonperturbative Ito's prescription is therefore to 
replace the fuU propagator Gil'^\t,t') by 

(0(i', x'mt, f)), = {4>{t' + e, x')^{t, x)) (43) 

with e — > 0^ , that is to multiply it in Fourier space by the regularization factor 
exp(— iecj'): 

(0(c.',g')0(^,g)), EE e-""'(0(L.',g')0(L.,g)). (44) 

Note that when this function is evaluated in a uniform field configuration, it becomes 
proportional to S(uj + lo') in which case exp(— iecj') = exp(zeci;). 

Let us finally prove the following theorem which is a consequence of Ito's 
prescription: if the bare action of a non-equilibrium system and the regulator 
term ASk are proportional to -0, then the effective average action r„['0,'0] is also 
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proportional to ■(/'• The proof is as follows. Let us rewrite Tj^ij),^] as the sum of a 
term independent of and of a term proportional to V'" 

T.liPM = ^ (^(VXx)) + Vi(x)r(o^i)[x;V',Vi]) • (45) 

Then: 

/ 14(V'(x))-r„[V,0]. (46) 

J X 

By hypothesis, Vk=a = since F^^a = S (see Appendix 8.1) and S is proportional to 
ip. Let us now suppose that remains zero down to the scale kq- Then at kq — dn^ 
Vkq-cLk is vanishing if and only if the contribution to V^o-dK coming from the flow 
between kq and kq — dn is zero. This contribution is calculated from 



^0 



d^l Kc(V'(x))u„ =a«r,[VAO 

= iTry a,i?„(q)G,[q,-q;V,Vi = 0] 



(47) 



In this equation, is computed by inverting the matrix (r^o +^ko)- As V^.^^i^) — 0, 
the element Tk^ [^.,^^\ of is vanishing at = 0. Moreover, recalling that, 
by assumption, AiS^ does not have any (jxp term (that is, as in definition (105) of 
Appendix 8.2, = = d^.R'^), the matrix element of the second row and second 
column of is vanishing. Therefore, one obtains: 

2 7q (43) 
(e-"G(i'i)[-q,q] + e^^"G(i/)[q, -q]) . 

Now, because of the causality of G^ '^^[— q, q], its poles lie in the upper half complex 
plane of w. The regulator term exp(— iew) in the integral of the first term allows for 
the closure of the integration contour by a semi-circle at infinity in the lower half plane 
without changing the value of the integral. From the residue theorem, we conclude 
that the integral over u! is vanishing since no pole is enclosed in the integration contour. 
The same holds true for the integral of the second term and we therefore conclude 
that Vug-dK is vanishing. By iteration, at all scales, 14 = 0. 

Note that, had we neglected the regulator terms exp(±iea;) in the integrals above, 
we would not have found a vanishing flow for and we would have concluded 
incorrectly that this term was generated by the RG flow. Note also that once F^ 
has been rewritten as in Eq. (45) with = 0, its flow is entirely determined by 
that of Fk°'^^ which is better behaved because it involves two propagators - as can 
be checked by taking the derivative of Eq. (37) w.r.t. ip. It turns out that in all the 
systems studied so far (e.g. [5, 6]), the resulting integrals were unambiguous and the 
regulator terms exp(±zea;) not necessary. We thus conjecture that the shortcut for 
Ito's prescription is to impose that F^ is proportional to ip (at least if F^ is an analytic 
functional of Let us now show how the NPRG formalism can be used to study 
Model A. 
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5. The derivative expansion at order two applied to Model A 

Model A describes the purely dissipative relaxation of a non-conserved scalar field 
(j){t,x) with Ising symmetry. It corresponds to Glauber dynamics (single spin flips). 
The model is defined by the Langevin equation (7) with N{(f)) = 1 and F deriving 
from the standard (equilibrium) Hamiltonian 

^[0] =^ J(V0)2+l/(0) with Vi4>)^^-4>^ + ^(l,\ (49) 

with the usual Z2 symmetry. The action of the model hence writes 

s[4>, i]= j^ X^ - v2<^ + y'(0)) - 0'} . (50) 

When approaching the continuous phase transition, the relaxation time of the order 
parameter starts diverging, which reflects the critical slowing down of the dynamics. 
Besides the static critical exponents v and r] of the Ising universality class, the critical 
dynamics is characterized by the dynamical exponent z which relates the divergences 
of the relaxation time r and of the correlation length ^ in the vicinity of the critical 
point as T ~ ~ |r — Tc\~^^ where Tc is the critical temperature. 

In the long-time limit this relaxation-towards-equilibrium model shows a time 
reversal symmetry that can be expressed as an invariance of the action (50) under the 
following field transformation [19, 20] 

{t -t 
t ^ t (51) 
4> 4>- dt(p. 

Indeed, on the one hand, the equilibrium part F = SH/dcf) of the action (50) is invariant 
on its own under the transformation (51) since the additional term cx dt4>5'H/5(l) 
vanishes upon time integration in the stationary regime. On the other hand, the 
time-evolution (j)dt4) and the noise (j)^ parts arc not invariant on their own but the 
combination 09f(/> — 4>^ is, since the terms generated by the transformation (51), that 
are proportional to {dt4>)'^ and 4'9t4'^ cancel out or combine to give back the original 
terms of the action. We refer the reader to Ref. [19] for a detailed and general study 
of the field-theoretic formulation of the time reversal symmetry. 

The computation of the static critical exponents of the Ising model has been 
performed in all dimensions using the derivative expansion [14, 5]. We now show 
how to compute the dynamical exponent z within this framework (using a richer 
approximation than in the former calculation of Ref. [ ]). 

The exact flow equations of the correlation and response functions do not form a 
closed set of equations since the flow of rlf-* involves rif-* and rl'^-' whose flow equations 
involve in turn rlf'' and tI^\ etc. As already mentioned, the derivative expansion is 
probably the simplest and most popular approximation - in particular for the study 
of critical properties - to close this hierarchy of equations. It amounts to proposing 
an Ansatz for under the form of a gradient and time-derivative expansion that 
corresponds to an expansion of all correlation functions in terms of their frequencies 
and momenta. 
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The Ansatz for T ^ at order one in time and two in space derivatives writes [20] 



X. 



and the initial conditions of the flow are 



(52) 

V' , Za ~ 1 ^ X\. Let us briefly justify 
this form: the time-reversal symmetry requires the Ansatz to be invariant under the 
transformation (51). It first implies that the ip'^ and ^pdtip terms renormalize in the 
same way. These terms henee bear in (52) the same coefficient Xi^^ip). It then implies 
that the term linear in -0 is invariant on its own, which imposes that it derives from 
a functional Hi^lip). We naturally choose for this functional the standard Ansatz for 
the equilibrium Ising model Hk = J ^k('0)(V?/')^/2 + [/^(V') [20]. Finally, note that in 
[20], only a field-independent running coefficient X^iip) = X^ was considered (leading 
order for the critical exponent z), whereas we here allow for a field dependence of 

x.ii^). 

The definitions of the functions involved in Eq. (52) are 



x.W 



= FT 



*=(V'.0) 



(57/>(x)(5V'(y) 

a.^FT ( r« 

(5-0(x)(5V'(y) 



i/=0,p=0 

* = (V':0) 



jy = 0,p=0 

u=0,p=0 



(53) 



where tp is a constant (in time and space) and where FT(.) means the Fourier 
transform. By translational invariance the momentum and frequency are vanishing 
in the first line and the second last two lines only depend on one momentum p and 
one frequency (see conventions of Eq. (97) - the trivial 2tt and 5{.) factors are not 
made explicit). In the spirit of the derivative expansion, the renormalization functions 
are computed at vanishing external momentum and frequency since within this 
approximation only the small momentum and frequency sector is correctly described. 

Their fiow follows from the flow of the one- and two-point functions, derived from 
Eq. (37). For f/;^ for instance, it is given by 



FT 













*=(V,o)) 










(54) 



where the 3-point function is put in the 2x2 matrix form F = dV)^ /dip and 

where a„ = d^Rj^d / dR,^. Taking the appropriate functional derivatives of (52) and 
evaluating the result at the uniform and stationary field configuration '^{t, x) = ("0, 0) 
one finds in Fourier space (see Appendix 8.2 for notations): 



-2X' 



(55) 



The propagator Gk in Eq. (54) is obtained by inverting the 2x2 matrix {t^k^ + R^^ 
evaluated in \['(i, x) = {ip, 0) (or equivalently from the general expression (113,114) of 
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Appendix 8.2). We find: 



V') 



1 



2X«(V') 
h{q) - iujXi^iiJj) 



h{q) 



- iLuX^{i}j) 




(56) 



with h{q) = Z^{il')q'^ + Rf^{q'^) + U'^{tp). More precisely, implementing the 
regularization advocated in Section 4, the off-diagonal terms of this propagator should 
be replaced by 



h{q)^ + {X,{^)ujf 



(57) 



After Fourier transforming Eq. (54) and inserting the expressions (55) and (56), the 
matrix product, trace and integral over the internal frequency ui arc straightforward. 
The resulting flow equation for C/^ is expressed in term of a single integral over the 
internal momentum q. The flows of Z„ and X„ are computed using Eq. (38), inserting 
the expressions for G^, F^'^^ and rif-' (see Appendix 8.2) and taking the appropriate 
derivatives with respect to the external momentum and frequency respectively as in 
(53). Once again, the internal frequency integral can be calculated analytically, such 
that the flow equations for and X^^ involve each only one remaining integral over 
the internal momentum q. 

Since we are interested in the scale invariant (fixed point) regime, we introduce 
dimensionlcss field and renormalization functions 



(58) 



where the running coefficients = Z,^{ipo) and = Xi^{ipo) arc defined at a fixed 
normalization point ipQ. In the critical regime, these running coefficients are expected 
to behave as power laws Z„ ^ n-vziK.) ^^^^^ ^ f^-vxM -^^ith ryz(K) — — k^kIuZk 
and similarly for %(k). The critical exponents 77 and z can be expressed in terms of 
the fixed point values of %(/*) and %(«;) a.s r] = r/* and z = 2 — 77* + 77*. 

As already mentioned, the three flow equations for U'^^, Z^ and X^, involve each a 
single integral over a momentum variable q. It is convenient to further introduce 
the dimensionlcss square internal momentum y = q^ j and the dimensionlcss 
cutoff function r(y) = {Z i^q^)~^ RK,{q^ / k^) . The flow equations for the dimensionless 
renormalization functions u', z and x are the sum of a dimensional part that comes 
from the change of variables (58) and of a dynamical part that comes from the 
integration of the rapid modes (previous calculations). The dimensional parts for 
u'jZ and x are respectively: 



dsU 
dsZ 
d.x 



dim 
dim 
dim 



= i (-(d + 2-?7,)u' + (-2 + d + 77z)^?i") 
= ?7x.T + i(-2 + d + ?7z)?Ax' 



(59) 
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and for the dynamical parts we find: 



dsX 



dyn 
dyn 

dyn 



= dyy^ - 



= Vd / dyy^ ^-^ 



/2 / ,.iyh'^ 2 



Vd / dyy^' 



s 



d h I d 
4/i2 



X 

h 2 



(60) 



with 9, = - 2'^7r'^/2r(d/2), /i(y,V) = y(z(V') + r(y)) + u"(V'), /(y,V) = 

yz'{tp) +u'"{ip) and s(y) = —Vz^iy) — 2yr'{y). As expected, x does not contribute to 
dsu' and 9s z which are the standard equiUbrium equations of the Ising model at second 
order in the derivative expansion [14, l ')]. The numerical study of these equations in 
dimensions two and three, which does not present any serious difficulty, will appear 
elsewhere. 



6. Supersymmetry 

It is well-known that the field theory associated with the Langevin equation of a model 
that relaxes towards thermodynamic equilibrium is endowed with a supersymmetry. 
Here we show that the supcrficld formalism that follows from this property leads 
to a rather simple NPRG formalism and, most importantly for us, is free from any 
ambiguity coming from the continuous time-limit. It is therefore a benchmark for 
testing our regularization prescription in Ito's discretization. 



6.1. Time reversibility and supersymmetry of the action 

We here follow the first route sketched out in Section 2.2 to treat the determinant 
(13) which consists in rewriting it using Grassmann variables 77, fj as in (14). We here 
consider that N{(j)) = 1. Once the Jacobian is exponentiated, one can finally integrate 
in (11) over the Gaussian noise distribution and deduce the generating functional of 
correlation and response functions 

Z[j,j] = J V(l)V[i4>]VT]Vfie-^^'''^^-'^-'^^+^''^^'^+''^ (61) 



with 



5[(A,</.,7?,^] = J ^^cj>{dt^ + niF{<P))~n2^''-v[dt + ni^^]fj}. (62) 
If F derives from a Hamiltonian: 



then the Langevin equation (7) (with iV = 1) corresponds to a dynamics that leads 
at long time to thermodynamic equilibrium. In this case, the noise strength is related 
to the relaxation rate by Einstein relation fli — fl2 (where ksT is set to unity), and 
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this coefficient can be convenientiy scaied away. Ttiis equiiibriuin property impiies 
that the action S possesses a supersymmetry, and it admits a compact form in the 
superspace in term of a (bosonic) superfield [12] 

X, 9, e) = $(n7) = (j){t, x) + ri{t, x)0 + 9fj{t, x) + 9e4>{t, x) (64) 

with Tu — (t, X, 9, 9) and where 9 and 9 are two anticommuting Grassmann variables 

{9,9} ^9^ = 0. (65) 

The integrals over these variables are defined as 

d9 = jd9 = Jd99 = Jd99 = l. (66) 

One introduces the differential operators 

D^d0-9dt, D = dg, with {D, D} = ~dt. (67) 
such that the generating functional takes the simple form 

Z[J]= j e-'5[*l+/» (68) 

where = J d'^x dt d9 d9, with the action 

^ I D^D^ + ($) (69) 

and the supersource J = j + + 7^ + 99]. The generators of the supersymmetry 
trnasformations are 

Q = de Q = dg + 9dt with {Q,Q}^dt. (70) 

One can check that the action (69) is invariant under the infinitesimal transformations 
(5$ = e(5$ and (5$ = eQ^. The operators D and D are hence covariant derivatives 
for the supersymmetry [12]. 

The correlation function between two superfields C{'cdi,tu2) = (<I'(roi)$(n72)) 
encodes all two-point correlation and response functions, and the supersymmetry 
encompasses all the dynamical symmetries. Indeed, the correlation function C then 
vanishes under the action of each of the three generators Q12, {Q, Q}i2 and Q12 (where 
Q12 = Qi + Q2 = dg-^ + dg^ and similarly for the others), which implies respectively 
causality, time-translational invariance and time-reversal symmetry, and yields in turn 
the fiuctuation-dissipation theorem. The supersymmetric formalism thus provides an 
elegant and powerful framework to treat equilibrium dynamics. 

Let us emphasize that reaction-diffusion systems do not, in general, show time- 
reversal symmetry. This is particularly spectacular for phase transitions to absorbing 
states. The directed percolation action (6) is invariant under the 'rapidity symmetry' 
(j){t) — >■ —(j)*{—t), (j)*{t) — 7> —(pi^—t) (which is a characteristic feature of the models 
belonging to this universality class), but it is not invariant under the transformation 
(51) and thus does not show time- reversal symmetry. For this reason, the field theory 
of this model cannot be recast into the superfield formalism with a simple action such 
as (69). 
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6.2. The NPRG in the supersymmetric formalism 

It is almost straightforward to render supersymmetric the NPRG formahsm. Here 
again a scale dependent generating functional ^^[J] is built from Z[J] defined in Eq. 
(68) by adding a quadratic cutoff term 

Z4J] = J p$e-5[*]-A5.[*]+/„ ./* (71) 

where $ is the supcrficld (64). The cutoff-term now writes 

A54$] = i/ <Piw)R,{w-m'Mw'), (72) 

where Rk{zu — w') = i?K(x — x')5(0 — 9')S{9 — 9') and where the Dirac function for a 
Grassmann variable is simply the identity: S(9) = 9. The functional is defined as 
the (modified) Legendre transform of log Z„ [J] : 

[*] + log [J] ^ I J-^- AS^ [*] . (73) 

with = (<I>). The exact flow equation for is formally identical to the one for a 
scalar field theory: 

a«r, = ^ / d^R^{w-w')G^\w,w';^\ (74) 

with Gk. = (^ri^' + and rl"'' the n-th functional derivative of with respect 

to the superfield '^{w) (see Appendix 8.3). The whole formalism is thus extremely 
close to the one for equilibrium theories, the only difference being the Grassmann 
dimensions that encompass the out-of-cquilibrium aspect of the model. Let us now 
show on the example of Model A how this formalism can be used in practice within 
the derivative expansion. 



6.3. Model A within the supersymmetric NPRG 

The supersymmetric action of Model A is given by Eqs. (69,49). As already stressed, 
the derivative expansion consists in constructing an Ansatz for that captures the 
low momentum and frequency sector of the model. In the supersymmetric formalism, 
it simply consists in a scries expansion of in powers of V and of the covariant 
derivatives D and D that is both supersymmetric and Z2-invariant. At order two in 
derivatives, it writes: 

Fk = y" I^X4-^)D^D^ + i Z,(*) (V*)2 + ;7k(*)| . (75) 

By integrating over the Grassmann variables in (75), one can check that the resulting 
Ansatz coincides with the non-supersymmetric Ansatz (52) with the same functions 
Uk, and X^. The flow equations for these functions are obtained by taking the 
appropriate functional derivatives of the flow equation (74). The effective potential 
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is defined from F^^^ evaluated in a uniform and stationary superfield configuration 



r 



(1) 



Hence the flow equation of (the first derivative of) the potential is given by 

5 



(5^'(777) 



(76) 



(77) 



Note that the cutoflt function is here chosen time- independent (see Section 4). As 

(2) 

usual, the full propagator in Eq. (77) is obtained by inverting (F^ + Rk) evaluated 
in the constant superfield configuration ^ = ijj (see Appendix 8.3 for more details). 
Wc find in Fourier space: 

G^{uj, q, Oi, 82) = ai(uj, qjS^iQi - 63) + a2(t^, 9)^^(61, 83) + 03(0;, q). 



with Q = 
and 



(78) 

52(61 - 62) = S{9i - 92)6{e, - 62), ^2(61, 82) = {Oi - e2){9i + 62) 



(79) 

C/^(V). Note 



aii(^,q) = h{q)/P{q,uj,ip) 
02(0;, q) = -iujXi^{'4>)/P{q,uj,'ip) 
a3iLU,q) = 2X«(V)/P(g,c.» 

where P(q,w, V) - h^io) + {X4i^)ujy and h{q) ^ Z^{i^)q^ ^ 
that within the supersymmetric formalism, the equivalent of the flow (47), i.e. Eq. 
(74) evaluated at "^{w) = ip, is vanishing by virtue of the Grassmann directions. 
No regularization is needed, the Grassmann directions automatically encode the 
appropriate causality properties. 

The Tk function evaluated in the same configuration is obtained by taking three 
functional derivatives of the Ansatz (75) and evaluating the result at ^ = ip. One 
finds in Fourier space, using the conventions of Eq. (97), (see also Appendix 8.3) : 

^('^(qi, 81, q2, 82, 83) = -2X',{i^) (52(61 - 62) + 52(61 - 83)) 

-x^(^) ((01 - 03)(^i - 02) + {Oi - e2){h - h)) 



+{u^^\i>) - z'MW? + 92' + 9i-<f2))<5'(ei - 82)52(61 - 63) 



(80) 



+2i(wi + w2)A:^(V')(^i - ^3)(ei - Q2){QxB2 + ^3^2) + 2 o 3. 

where, by translational invariancc, ^3 = — oji — l>J2. Once Eq. (77) is Fourier 
transformed in space and time, it is straightforward to insert in this equation the 
expressions (78), (79) and (80) and to perform the integrals over the Grassmann 
variables, which cancels out most of the terms. The non- vanishing contributions yield 
the flow of [/^ (?/'), in which the integral over the internal frequency a; can be achieved 
to express this flow in term of the only integral over the internal momentum q that 
remains. The result for the dimensionlcss potential is identical to the one obtained in 
the nonsuperymmctric formalism, Eq. (60). 

(2) 

The flow equations for Z„ and can be deflned from the flow equation of F^ 
evaluated at the constant fleld conflguration ^'[nj] = 



ap-2FT (d^Tl^\mi,W2: 
a,.FT (d^T^^\nJi,nj2) 



(81) 
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with the same conventions as in Eq. (53) and with the flow of given by 

5.r(2) = d.R. (-ir(4) + T^^G.T^^ G., (82) 

where the arguments ({tUi}; ■(/') of the vertex functions and of the propagator are 
imphcit. The expression for Fj; is derived in Appendix 8.3. Once again, after 
inserting (78), (79) and the expression for rlf'' into Eq. (81), after performing the 
integrations over the Grassmann variables and over the internal frequency, the flow 
equations depend on a single momentum integral and the equations thus obtained are 
identical to those derived in the nonsupersymmetric case, Eqs. (60). This constitutes 
a nontrivial check of the consistency of the regularization of the propagator proposed 
with Ito's discretization. 

7. Conclusion 

This paper was devoted to presenting in detail the NPRG formalism for non- 
equilibrium stochastic models in statistical physics. Special emphasis has been 
dedicated to the inspection of the consequences of Ito's choice of discretization within 
this formalism. The main objectives were three-fold i) provide a reference including all 
the specificities of NPRG for non-equilibrium systems ii) set up the NPRG framework 
to deal with supersymmetry, iii) identify the relevant prescriptions in NPRG related 
to Ito's discretization and device a systematic way to enforce these prescriptions. 

Regarding ii) we showed in Section 6 that the supersymmetric formalism can 
be simply incorporated within the NPRG framework, and we derived in this context 
the NPRG flow equations for Model A at second order in the derivative expansion. 
(Of course, this supersymmetric version of NPRG is mostly useful when time-reversal 
symmetry holds.) As for iii), going back to the discrete version of the field theory 
ensueing from either a Janssen-de Dominicis or a Doi-Peliti transformation, we 
analyzed the causality properties of the response functions and showed that these 
properties are preserved under the NPRG flows. We proposed a simple prescription 
- a regularization factor e^"^" in the propagator - to enforce automatically the 
consequences of Ito's prescription within the NPRG flows |. We illustrated this 
procedure on the non-supersymmetric calculation of the NPRG flow equations for 
Model A, which identify with the ones obtained in the supersymmetric version and 
thus constitute a non-trivial check of the consistency of the proposed prescription. We 
emphasize that, of course, the whole non-supersymmetric framework presented here 
also applies to genuinely out-of-equilibrium problems. 

To conclude, let us emphasize once more that, as the continuous-time field theories 
for non-equilibrium models (derived either from a Janssen-de Dominicis or Doi-Peliti 
formalism) are ill-defined on their own, it is crucial to be aware of the implicit dis- 
cretization chosen and of its implications. This applies in particular within the NPRG 
framework. This is the first time that an exhaustive account of the specificities of the 

X Note that this prescription is very close in spirit to what is encountered within the second- 
quantization formaUsm, for instance in sohd state physics, in the presence of a first order time- 
derivative. Once the Hamihonian is put in a normal ordered form, partition functions can be written 
in term of a functional integral expressed in the coherent states basis. When taking the continuous- 
time limit, the ambiguity of the value of 9{t = 0) arises, as in the Doi-Peliti formalism, and is solved 
in a way similar to the one advocated here, via a regularization factor exp(±iei.ij) in the propagator. 
See, e.g. G.D. Mahan, Many- Particle Physics, 2nd. edition, (Plenum, New York, 1990). 



General framework of NPRG for NESS 



22 



NPRG methods for non-equilibrium systems is provided, including a detailed anal- 
ysis of discretization problems in this context. This contribution, together with the 
regularization proposed, will serve as a reference for future NPRG calculations in non- 
equilibrium models. 
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8. Appendix 

8.1. The limit k — > A.' a proof that Fa = S 

We prove in the following that if the cutoff function diverges at k = A, the effective 
average action Pa becomes the bare action at this scale. For simplicity, we consider in 
the following only the case of a frequency-independent cutoff function i?^ . The proof 
is rather different for the field theories derived from a Langcvin equation (Janssen-de 
Dominicis formalism) and from a master equation (Doi-Peliti formalism). We thus 
study both cases separately. 

8.1.1. Master equations 

The Doi-Peliti formalism for master equations naturally leads to field theories 
written in terms of two fields </> and 4>* that are complex conjugate. It is convenient 
to write 

(j) = (j)i+ i(t>2 (83) 
with 4>i real and to couple these fields to real sources Ji in the following way: 

Note that with this choice the fields (f> and 4>* are coupled to independent and real 
sources: 2{Ji(j)i + J2i(j)2) = {Ji + J2)4' + (^^i ~ J2)4'* ■ From the definition (84), it is 
straightforward to show that Z* ~ which is thus real. Note that this property 
follows from the fact that, although complex, the action 5 is a polynomial in 0i and 
i(/)2 with real coefficients (the reaction rates). It is also straightforward to show that 

= and (i<^2)* = («02)- (85) 

We define ipi = V'2 = (*02) and — {(f)), i/j — {(j)*) that are all real and 

independent: ip = ipi + ip2 and ijj — ipi ~ 'ip2- We choose in the following the natural 
and convenient cutoff term A5k[(/', (/<*] = AiSk[(/)i, i(/)2] 

AS, = \j^R.{q^)MW{~^) - ^^i?«('f')('^i(q)0i(-q) + </'2(q)02(-q)) (86) 
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which we moreover take frequency- independent. We symbohcahy note it 1/2 j Rf^{(j)\ + 
</)|). We then define the effective average action as 

T^[il^i,^2]+W^[JuJ2] = j {Jii^i + J2^2)-\j ^^'(x)i?«(x-x')^(x'). (87) 

with Wk. = logZ^. Note that this definition leads to a minus sign in the last term 
when written in terms of ^/^i and "02: 1/2 / Ruii'i ^ ^l)- It is now easy to prove from 
Eq. (87) the following relation: 

g-r„[)/'i,i/'2] ^ J 'p^^qy^^^-S[cl>i+i>i,i4>2+i'2]-AS^[^iA4>2] + J ^<t>i + ^i<t>2 ^ ^gg^ 

The cutoff term A5„[0i, i(/)2] in Eq. (86) is positive and formally identical to the 
equilibrium one for a 0{2) model. This implies that, as usual, if R^ becomes 
infinite when k — > A the only configuration of 0i and 02 that contributes to the 
functional integral in Eq. (88) is 0i = 02 = 0. This leads to the expected result 
rA['0ij'02] = = ipi,i(f>2 = ip2]- Note that since the action 5 of a reaction- 

diffusion system is proportional to 0*, this result shows that Fa is proportional to 



8.1.2. Langevin equations 

The situation is different for Langevin equations because the fields and are 
respectively real and pure imaginary (note that in Eq. (11) the integration over is 
performed on the imaginary axis). In Ito's discretization where Grassmann fields are 
not necessary, the action corresponding to a generic Langevin equation is quadratic 
in (once the average over the noise distribution has been performed). It writes: 

S = J 4>K{cb)~$'N^{cf>), (89) 

where i^(0) = 9*0 — F{4)). The generating functional of correlation and response 
function writes 

Z^[j3] = j V(t)Vi4>e-^^'f''^^-'^^''^'t''^^+^ (90) 

with both j and j real. This expression is formally very close to Eq. (84) up to the 
fact that the cutoff term writes 

AS, = I i?«(g2)0(q)0(-q) (91) 
and is thus also pure imaginary. The definition of the effective average action is 

r^V, ^] + W,[J, J]= f (j^j + jv3) - / ^(x)i?4x - x')V'(x') (92) 

Jx V x,x' 

where "0 = (0) and ip = (0). Once again, upon rewritting as i0 with real, 
it is straightforward to show that Zn[j,j], ip and ip are real. Following the same 
steps as above one can derive the same relation as Eq. (88) with 0i and i(j>2 replaced 
respectively by and 0, and 4'i a-nd "02 replaced by ip and ip- Now, since S is quadratic 
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in it is possible to perform the functional integral over this field (taking into account 
the fact that it is imaginary). We obtain: 

X exp (^J ^(f> - J ipK{(j) + V) - 4'^N\(j) + V-)^ . 

(93) 

Once again, if diverges when k — >■ A only the configuration = contributes to 
the functional integral above. By expanding Fa as a power series in tjj: 

Ta = / Aoi^) + MiW + i'^A^ii^) + ... (94) 



we find from Eq. (93) evaluated at </) = 0: Aa{^) = 0, Ai(?A) = N^{^P) , A2{^j) = G{4') 
and all An{ip) = for n > 2. We thus conclude that rA[-0,-0] = S[4i ^ 4\4> = V'] as 
expected. Note that this result proves that rA['0,'0] oc ip- 

8.2. Formulas for the NPRG propagator and vertex functions 

In the nonsupersymmetric case, we only need space and time coordinates and we 
define x = (<, x) and q = (lo, q). Our convention for the Fourier transform is 

/(x) = / /(q)e'(^~-^~-'^*) (95) 

where the same symbol is (abusively) used for the function and its Fourier transform 
and where 

For a translational invariant function /(x, y) which depends only on the difference 
(x — y), we define: 

/(q, q') = (27r)'^+l5'^+l(q + q')/(q), (97) 

without again changing symbols for the function. The generalization to functions of 
n momenta is straightforward: /(qi,...,q„) = (27r)''+M''+^(X;i qO/(qi> • ■ • 7 Qn-i)- 

The matrix of two-point connected correlation (and response) functions is 
by definition: 

W^^^ fx x'- J] - (98) 

(1) 

The same definitions hold for the vertex functions F)..^j[x,x'; by replacing Wk, by 
Fk and J by ^ . From now on, the k and J (resp. 'J') dependences of W (resp. F) are 
not systematically indicated to alleviate the notations. We introduce the alternative 
notation, e.g. 
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where the first (resp. second) exponent in (1, 1) refers to the number of derivatives with 
respeet to j (resp. j). The relationship between the two definitions of the connected 
two-point functions is therefore: 

(100) 
that is 

W^(2'^[x,x'] = W^(i'i)[x,x'] = W^if)[x',x] 

W^i^i^[x,x'] = M^(2,o)[x^x/] ^ VFi2^[x,x'] = iy("'2)[x,x']. ^ ' 

The same kind of relations exist between the two-point vertex functions F*^"'") (x, x'), 
n + h = 2. In Fourier space, this imphes: 



\ ^[q',q; J] 1q,q'; J] 



(102) 



When these functions are evaluated in a uniform field configuration, this matrix 
becomes proportional to (5'^"'"^(q -I- q') and using the convention of Eq. (97) yields: 

M/(2)(q;j)^ f for J uniform (103) 

and similar relations for Fk (q; ^) when ^ is uniform. 

The modified Legendre transform F^ (the "efi'ective average action") is defined 

by 

r.[*] + W,[J] = /" * J.$ - i /" ^ **(x).i?,(x - x').*(x). (104) 

■/X •/ x,x 

In Fourier space: 

"-(o) -(«;•(-<,) SS) 

where we have taken i?^"(q) = 0. The true Legendre transform of is therefore 
Fk[*] + ^ /**.-Rk-*- We thus find that the full field-dependent propagator is: 

G„ = (f(2)+i?,)-i = M/(2) (106) 

since 

/ W^i2)[x,y;J].(fl2)+i?,)[y,z;vI/]=5'i+i(x-z)l2 (107) 
where I2 is the 2x2 unit matrix. In Fourier space this relation writes: 

/ VF(2)[p^q. (r(^2) ^^^)[_q^p,.^] ^ (2nf+^5''+\p + p')h. (108) 

The exact flow equation for F^ writes in Fourier space: 

d.^^m = \^^ j 5Ki?4q)-G«[q, -q; *] (109) 
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with defined by Eqs. (102,106,107). Note that this equation can be conveniently 
rewritten as 

d^T^m = i 9, Tr^ In (f^^^) + ^ (i^g) 

where = di^R^d / Or^, which is frequently used. 

It is in general sufficient to consider the flows of the vertex functions evaluated in 
a uniform field configuration (for the calculations performed within the derivative 
expansion for instance). In this case, the equations simplify since the problem becomes 
translational invariant. The flow of the two-point functions in the configuration '^u 
writes: 

9.ri%(p;vI/„)=Tr j a.i?(q).G(q). 

rf ^ (p, q, -P - q).G(p + q).f f ' (-p, p + q, -q) - ^ f If (p, -p, q, -q)) .G'(q) 

(111) 

where, in the uniform fleld configuration and using the convention of Eq. (97) : 

G«(q, *„).(f (2)(q, *„) + i?,(q)) ^ (112) 
and thus, omitting to write the ^„ dependences: 



r 1 / -rP'(q)-i?!!^(q) r^^'(q) + i?ii(q) 1 ^i^^ 



with 

P(q)= -r^°)(q)(r^^)(q) + i?02(q)) 

+(rl'^' Vq) + i?iH-q))(rL'^'^(q) + R]H^))- 

The three-point functions appearing in Eq. (Ill) are defined by 



(114) 



r. [^.y.^;^= (115) 

and in (111), they are evaluated in a uniform field configuration and Fourier 
transformed. The same kind of definition holds of course for the four-point functions. 
In a uniform field configuration, the precise relationship between the s and the 



F^^Vn n n - f ^1^'"^ (P' q' " q) tI'''^ (p, q, -p - q) , 

Tk (p, -p - q, q) F^ ' ' (p, q, -p - q) 



and 



f(3)^„ „ „_„N ( r- (q'-p-q.p) r«- ^(q,p,-p-q) , 

r!. n-p-q,P,q) Fi, '(p,q,-p-q) 



and for the four-point functions: 



f^^Vr, r. n - I ^'^ HP, q, -q, -p) HP, q, -P, -q) 1 /noN 

1 12 lP'-P>q:-q} - j.(2,2), . -p(1.3)/ X (11^'' 

r;. (p, -q, -P, q F^ ' ' (p, -p, q, -q 
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f^^Vn nn - ( ^^''^ (p, -p, q, -q) \ 

J- 11 IP, -p,q, -qj - X -p(2,2). X U^yj 

V ' '(p,-p,-q,q) rj, ' '(p,-p,q,-q) ) 

fWfn _n o _o^ _ rp)(q'-q'P'-p) rL''^'(q, p, -p, -q) \ /,2n^ 
\ (-q,p,-p,q) (p,-p,q, -q) / 

8. 3. Formulas for the supersymmetric formalism 

In this appendix the expressions of the supersymmetric propagator GkI'cui, tu2', 5*] and 
vertex functions ri?' and T^^"^ are derived. For the propagator, one needs to determine 
the inverse of (vl^^ + Rk^ . Taking two successive functional derivatives of defined 

by the Ansatz (75) with respect to "^{wi) and '^{w2) and evaluating the resulting 
expression at uniform and stationary superfield ^'(tu) = yields 



(fL'^ + i?,) {wi,W2) = [X,{-^){DD - DD)^, - ^«(^) 

+ R^{xi-X2) + U'^{^)Y'+^{mi-m2) 



(121) 



with (S'*+3(roi - W2) = S{ti - t2)S'^ixi - X2)S{0i - e2)S{9i ~ 62) and recalling that for 
Grassmann variables S{9) = 9. It is convenient to Fourier transform this expression 
in space and time only (and not in the Grassmann directions since there is no 
translational invariance in these directions). Using the convention of Eq. (97) - that 
is without writing the (27r) and 6{.) factors resulting from translational invariance in 
space and time - we find: 

(r^'^+i?,) (qi, 91,62) - [{Z^{i')q,^+R^{q,^) + U':W)S\ei~Q2) 

+iLUiX,{^b)S\ei,Q2) - 2X,(^)} 

(122) 

with e = {0j) and ^2(81-62) = s{ei-02)S{0i-e2), P{ei,&2) = (ei-^2)(^i+^2). 

Let us now determine its inverse G^. As ^F^f^ + i?^^ is diagonal in frequencies and 
momenta, so is G^. The inversion relation writes 

JdQs G«(c^, 9,81,83) (fI^) +7?^) (w, 9,83, 82) =-52(81 - 82). (123) 
We search for the inverse G^ of (tI^^ + Rk] assuming the same structure in the 



Grassmann directions, that is 

G«(l^, q, 81, 82) = ai{io, q)6^{ei - 82) + 02(10, q)PiQu 82) + 03(10, q). (124) 
The coefficients Oi can then be calculated using Eq. (123). One finds 

ai(uJ,q) = h(q)/P(q,u,il;) 

a2(uj,q) = ~tOjX^(i;)/P(q,uj,^) (125) 
a3(u;,q) = 2X^(^) / P(q,OJ, ^) 

where P(q,UJ,^) = h^q) + (X4^)ujy and h(q) = Z4^)q' + R^(q^) + U'^(^). 
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Let us eventually give the expressions of the 3- and 4- point vertex functions. By 
taking an additional functional derivative of the Ansatz (75) and evaluating the result 
at the uniform and stationary field configuration, one finds 

- Z',{^) [s-'+^w, ~ W2)Vl^5''+\w, ~ W3) 

+ 5'^+3(tni - tJ73)V|/'^+3(tUi - W2) 

+ Vj,,5'^+3(n7i - tJ72)V5,,5'^+3(tJ7i - n73) 

+ X'^[ij) [(5'^+3(tui - UJ2){DD - DD)^,5'^+^[vji - w^) 

+ (5'^+3(tui - VJ3)[DD - DD)^y+^{wi - W2) 
+ D„,(5'^+3(tui - TU2)5^i(5'^+3(n7i - 073) 

+ D^,5'^+3(wi - W3)D^,5'^+'^{wi - W2) 



(126) 



which in Fourier space writes 



rl')(qi,ei,q2,e2,e3) 



(c/f '(V) - Z'^mq,^ + qi + 91.92)) <52(ei - Q2)5\Ql 

3)) 



- 2X',{^l^){6HQ,-Q2) + 5\Qi-Q3)) 

+ x'^{-,p) {{Oi - 03)(^i - 02) + [91 - e2){ei 



+ 2i(wi + a;2)^;;(V')(^i - ^3)(^i - e2){h02 + hh) + 2 o 3. 

(127) 

where, by translational invariance, W3 = — cji — u)2- The expression of the 4-point 
vertex function follows from the same procedure, which yields in Fourier space 



63) 



rl''^(qi,ei,q2,e2,q3,e3,e4) 



= I [[/^''^ + z'^ (VO {ql + + (7i + 91 .92 + gi .93 + 92 .ga) 

X ,52(91 -e2)(52(ei- 83)52(0^ _ 04) 

- 2X;;'(V') [52(61 - 62)52(61 -63)+3f^4 + 2f^3 

+ X'^{^:) [52(61 - 62)(f?i - 04)(^i - ^3) + 3 o 4 

+ 52(6l-63)(^^l-04)(^l-^2) + 2o4 

+ 52(6i-64)(ei-02)(^i-^3) + 2f^3 

- x'^ii,) + c^2 + ^3)(^?i - 02){ei - - 64) 

X (^i^2^3 + ^2^3^4) + 3 ^ 4 + 2 O 4U , 



(128) 



where again, by translational invariance, 0-14 



-Ui - LO2- W3- 
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